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Oon the pictures before any combination is achieved. Some properties of said decomposition, in particular in the Fourier frequency 
domain, enable substantial reduction of the amount of computations. 
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(57) Abrege : Un calcul de reconstruction d'images en tomographic implique la decomposition de jeux bidimensionnels de pro- 
jections de cette image par une mSthode de decomposition en ondelette, dans laquelle le jeu de projections est separe* en imagettes 
exprimant respectivement r ensemble et les details a plus haute frequence des projections. La reconstruction s'effectue separement 
sur les imagettes avant qu'une combinaison ne soit faite. Certaines propri6t6s de cette decomposition, notamment dans le domaine 
frequentiel de Fourier, permettent de r&iuire sensiblement la quantity de calcul s. 
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PROCEDE DE RECONSTRUCTION D'UNE IMAGE A PARTIR D'UN JEU 
DE PROJECTIONS PAR APPLICATION D'UNE TRANSFORMEE EN 

ONDELETTE 



5 DESCRIPTION 

Le sujet de cette invention est un proc^de 
de reconstruction d'une image a partir de jeux de 
projections de cette image et en appliquant une 
transformation en ondelette . 

10 Les proced^s de tomographie consistent §l 

examiner un objet inanime ou un etre vivant par un 
reseau de detecteurs mobiles qui prend une serie de 
vues en tournant autour de lui . Ces vues sont des 
projections de la propriety permettant d'exprimer 

15 1' image (normalement une absorption d'un rayonnement 
traversant 1' objet ou une scintillation d'un corps 
radioactif inger£ par l'objet), c'est-a-dire des sommes 
de la propriSte le long de lignes traversant 1' objet et 
d6finies par la collimation des detecteurs. Chaque 

20 d^tecteur mesure une projection de 1 1 image k chaque 
vue. Quand un nombre suffisant de vues et de 
projections a §te pris, on precede Sl l 1 inversion des 
resultats pour obtenir la valeur de la propriety a 
chaque point de 1' objet ; cette inversion est 

25 comparable £ 1' inversion d'un systeme d' Equations de 
grandes dimensions et peut s'effectuer ouvertement par 
des mithodes algebriques, ou, plus frSquemment, par des 
m£thodes analytiques par lesquelles on applique des 
operations numeriques successives aux projections sans 

30 inverser directement le systdme. Un grand nombre de 
m^thodes existent, parmi lesquelles on citera celle qui 
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est d§crite dans le brevet frangais 2 615 619 qui est 
le premier de cette 6quipe de recherche, et le brevet 
frangais 2 810 141 plus recent qui prSsente quelques 
ressetnblances avec le procede qui sera decrit ici. On 
5 peut aussi citer un article de Grass paru dans Phy. 
Med. Biol, vol.45, p. 329, fgvrier 2000. Ces operations 
aboutissent k ce qu'on appelle la r£troprojection des 
mesures, c'est-a-dire au calcul de la valeur de la 
propriety prise k chacun des points de la ligne de la 
10 projection. 

II peut etre avantageux de travailler avec 
des r6sultats exprim^s dans le domaine friquentiel de 
Fourier, comme le prouve le deuxidme document. Des 
transformees numeriques d' autre nature ont aussi 6t6 
15 employees . 

On se reportera aux figures 6, 7 et 8 pour 
une description concrete et sch€matique d'un proc£d6 de 
prise de mesures tomographiques . Une source F de 
rayonnement et un systdme de detection 3 sont mobiles 

20 le long d'un cadre annulaire 2 £ des positions 
diametralement oppos6es, et un rayonnement conique 
originaire de la source F atteint le systeme de 
detection 3 apres avoir traverse I'objet 1 a 6tudier. 
La pi£ce essentielle du systdme de detection 3 est un 

25 reseau 4 bidimensionnel de detecteurs 5. Des 
projections R de 1 1 image tridimensionnelle de l'objet 1 
sont mesur^es par ceux des detecteurs 5 qui sont inclus 
dans un pourtour 15 de 1 1 11 ombre" de l'objet 1. Un 
nombre important de vues de ce genre est pris a autant 

30 d' angles 6 differents d 1 orientation du reseau 4 des 
detecteurs 5. 11 est frequemment considers, 
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conventionnellement pour les calculs, un reseau fictif 
4' de detecteurs 5 1 sur un plan de detection Pdet qui 
est parallele au reseau 4 reel et passe par le centre O 
du cadre 2. Des coordonnSes p et q sont definies pour 
5 identifier les detecteurs 13 d'apres leurs lignes et 
leurs colonnes. Des calculs de rearrangement, courants 
dans l*art, permettent de transposer les mesures d'un 
reseau 4 quelconque au reseau fictif 4' et d'appliquer 
les algorithmes de reconstruction k celui-ci. 

10 Le brevet fran<?ais 2 615 619 sera rappeie 

ici puisqu'il expliquait en detail quelles operations 
numeriques, notamment celles qui sont dites de filtrage 
et de retroprojection, permettent d'obtenir 1 1 image de 
I'objet k partir de ses projections ; des mSthodes 

15 differentes existent cependant. 

Le probldme auquel on s 1 attaque ici est la 
reduction du temps ou du volume des calculs pour 
inverser le systdme des mesures. 

On sait qu'il s'agit d'une des limitations 

20 les plus graves dans les proc6d£s de tomographie, et de 
nombreux proc6d£s nouveaux ont ete con?us pour le meme 
objectif que 1' invention, dont celui du deuxidme brevet 
cit6 • 

L'idee exposee ici est d 1 exploiter les 
25 proprietes particulidres d'une transformee numerique, 
dite decomposition en ondelette, des projections pour 
prevoir des emplacements de projections negligeables ou 
insignif iantes et de ne pas appliquer les calculs 
d 1 inversion a ces emplacements. Des perf ectionnements 
30 off rent encore une plus grande reduction des calculs. 
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Sous sa forme la plus g€n6rale, 1 1 invention concerne un 
proceed de reconstruction d'une image k partir de jeux 
de projections de cette image, comprenant 
successivement : 
5 - une serie de decompositions successives 

des jeux de projections en ondelette dormant des 
imagettes des jeux des projections, comprenant des 
images d'une approximation (AA) et des series 
successives (Dd, Dh, Dv) de details homologues de 
10 chaque jeu, 

- dans chacune des series et successivement 
pour des imagettes de details ayant un nombre croissant 
de points, une recherche de portions insignif iantes 
jug^es pour depourvues de contenu, et une recherche de 

15 portions homothetiques des portions insignif iantes dans 
chacune des imagettes de details qui suivent dans la 
serie ; 

- des retroprojections des imagettes des 
jeux de protection en imagettes de 1 1 image & 

2 0 reconstruire, avec une omission des retroprojections 
pour toutes les portions insignif iantes et toutes les 
portions homothetiques. 

et une combinaison des imagettes k 
reconstruire par une inversion de decomposition donnant 

2 5 ladite image. 

L'ordre des etapes (surtout les 

retroprojections precedant la recomposition) est 
essentiel pour obtenir les avantages de l 1 invention. 

L'art anterieur comprend un exemple de 

3 0 reconstruction d'une image par une decomposition en 

ondelette (US - 5 953 388 - A) , qui cependant est 
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appliquee Ik pour reconstruire une portion seulement de 
l 1 image, en exploitant la proprigte de "localite" de la 
decomposition, qui n'est guere sensible aux autres 
portions de 1' image, qu'on peut ainsi negliger dans les 
5 calculs. Le brevet US-5 841 890-A traite d'un sujet 
analogue et considere divers aspects de la 
decomposition en ondelette appliquee k la tomographie. 
La decomposition ne sert pas a acc616rer des calculs de 
retroprojection pour obtenir une image d 1 ensemble de 

10 fa<?on plus econome en calculs. 

On doit enfin mentionner le document US-6 
148 110-A qui decrit un procede a masques numeriques 
analogue a celui de l 1 invention, mais seulement pour 
comprimer un signal d' image sans calcul de 

15 retroprojection. 

L 1 invention sera maintenant decrite plus 
concretement et completement en liaison aux figures 
suivantes : 

la figure 1 illustre une decomposition 
20 d 1 image en ondelette, 

- les figures 2, 3 et 4 illustrent certains 
aspects de l 1 invention, 

la figure 5 illustre un organigramme 
resumant une realisation complete du procede, 
25 - la figure 6 illustre un dispositif de 

principe du proc6d6, 

- et les figures 7 et 8 illustrent plus en 
detail le procedS et certaines rotations employees. 

Nous commencerons par exposer la 
3 0 transformation d'un signal par une decomposition en 
ondelette. Plusieurs modeles d'ondelettes existent, qui 
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ont en commun d'etre comparables a un filtre passe-bas. 
Le signal est sipare en deux portions dont l'une, 
associee aux basses frequences, peut etre tenue pour 
une approximation du signal, alors que 1' autre, 
5 associSe aux hautes frequences, exprime plutot ses 
details. Une propriety des ondelettes est que les 
portions peuvent contenir chacun une moitie des points 
du signal si bien qu'il n'y a pas de perte 
d 1 information par cette decomposition. La decomposition 
10 peut §tre faite dans le domaine direct d 1 expression du 
signal ou dans le domaine de Fourier. 

Dans le cas de projections d'un objet 
examine sur un reseau generalement bidimensionnel de 
detecteurs, les projections peuvent St re groupies en 
15 jeux bidimensionnel s d'apres deux de leur coordonnles 
(en general p et q sur les axes d'un reseau de 
detecteurs). Toutefois, comme la figure 8 l'illustre, 
on considere le plus souvent des jeux de projections 
rearranges sur un reseau fictif de detecteurs. Dans 
20 I'exemple de la figure 8, on regroupe des projections 
Rx originaires d'un certain nombre de positions Fx 
successives de la source F du rayonnement, de maniere 
que les projections about issant a une meme colonne (k q 
constant) de detecteurs 5' du reseau fictif 4 1 soient 
25 originaires d'une meme position Fx, et aussi que les 
plans de projections soient tous paralleles jusqu'aux 
positions Fx : on a alors transform^ le probleme d'une 
geom^trie conique du rayonnement & une ggom6trie en 
eventails paralldles dont la resolution et plus simple. 
30 De plus, le reseau fictif 4' passe par le centre 0 de 
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rotation et appartient done au plan de detection Pdet, 
ce qui facilite aussi les calculs. 

L 1 invention pourrait encore etre appliqu€e 
a des rencont ructions de coupes de l'objet 1 a l'aide 
5 d'un reseau monodimensionnel de detecteurs (tous & la 
meme coordonnee p) . On pourrait encore appliquer le 
principe exposg ci-dessus du rearrangement des 
projections d'une geomStrie en event ail plan a 
geometrie parallele. 

10 Le traitement des mesures rearrangees se 

fait en suivant les lignes de detecteurs 5 du reseau 
fictif 4', success ivement pour tous les points. 

Un signal peut §tre decompose 

successivement en ondelette pour dormer plusieurs 

15 niveaux de resultats. Les nouvelles decompositions ne 
concernent que la portion de 1' ondelette qui donnait 
1' approximation du signal, la portion ou les portions 
qui donnaient les details etant pr6serv6es. 

Prenons comme exemple d'objet d'une 

20 decomposition en ondelette une image form£e de cinq 
cercles dont un cercle exterieur et quatre cercles de 
diametres differents, tous inscrits dans le premier. 
L' image de depart comprenait n x n points, et si on 
applique deux fois une decomposition en ondelette de 

25 cette image selon le principe precedent, le resultat 
est donne & la figure 1. 

La decomposition en ondelettes de 1' image 
donne un ensemble d'imagettes, dont trois sont plus 
importantes que les autres, comprenant n/2 x n/2 points 

30 chacune, et correspondent aux details horizontaux, aux 
details diagonaux et aux details verticaux de 1' image 
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de depart a grande ichelle ; elles sont notees 
respect ivement Dhl, Ddl et Dvl . Les details horizontaux 
de 1 ' image sont obtenus a partir des projections 
d' angles 0 proche de o ou n, les details verticaux k 
5 partir des projections d 1 angles 8 proche de n/ 2 ou 3rc/ 2 , 
et les details diagonaux a partir des projections 
d 1 angles interm€diaires avec les conventions de la 
figure 1. Le reste de 1' image est compose de quatre 
imagettes comprenant n/4 x n/4 points chacune et dont 

10 trois sont des imagettes de details horizontaux, 
diagonaux et verticaux k plus petite echelle, notees 
Dh2, Dd2 et Dv2, alors que la derniere imagette est une 
approximation notle AA de 1' image de depart. Si on 
appelle <|> et \|f les fonctions de decomposition d'une 

15 image ou d'une imagette en ondelette, la fonction <|> 
donnant 1 ' approximation et la fonction \\f les details, 
les fonctions a appliquer k 1' image de depart pour 
obtenir la decomposition de la figure 1 sont donnees 
par le tableau I . 

20 AA : * (x a ) <|> (x 2 ) $ (^) $ (y) 

Dv2 : <J> (x x ) <|> (x 2 ) <|) (y) (^), 
Dd2 : <|> ( Xl ) (|> (x 2 ) (^) *F (^), 
Dh2 : (|> (Xx) <|> (x 2 ) V (^) <J) (-y) 

DVI : <|> ¥ (X 2 ) 

25 Ddl : ¥ (xx) ¥ (x 2 ) 

Dhl : ¥ (xx) (J) (x 2 ) 

L 1 invention consiste effectuer la 

r€troprojection sur chacune des imagettes des jeux de 
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projections d6compos§es en ondelette et k combiner les 
imagettes retroprojetees par une inversion de la 
disposition en ondelette pour obtenir 1' image 
recherchee. La decomposition en ondelettes est propice 
5 a diverses simplifications qui acc61erent beaucoup les 
calculs. Ces simplifications sont faites entre la 
decomposition et la combinaison. 

La premiere d' entre elles concerne des 
filiations qu'on peut etablir entre des details 

10 homologues a des echelles differentes. On considere 
pour cela des series d' imagettes dormant des details de 
m§me nature. La figure 2 (qui illustre une 
decomposition d'une image ressemblant k celle de la 
figure 1, mais a trois niveaux de decomposition) 

15 illustre, pour les trois imagettes des details 
horizontaux Dhl, Dh2 et Dh3 , des portions homologues 
Jl, J2 et J3 qui occupent la m§me position et la mime 
superficie relative sur chacune de ces imagettes et se 
deduisent ainsi l'une de 1' autre par une homothetie 

20 geometrique dans leurs imagettes. 

On peut faire l'hypothese que, pour la 
plupart des images rencontrees en pratique (a 
1' exception notamment des images texturees) , si une 
portion telle que J3 a un contenu insignif iant , 

25 c'est-a-dire qui ne revele rien par rapport aux valeurs 
signif icatives de l'imagette, les portions homologues a 
plus grande echelle, ici J2 et Jl, seront elles aussi 
insignif iantes . 

Conformement & 1' invention, on commence 

30 done, pour les imagettes des details, k retroprojeter 
les details £ 1' echelle la plus petite, puis les 



FEUILLE DE REMPLACEMENT (REGLE 26) 



WO 2004/036288 



PCT/FR2003/050089 



10 

details aux echelles de plus en plus grandes. Un seuil 
numerique est appliquS aux coefficients de 1'ondelette, 
c'est-&-dire aux valeurs que prend la transformee dans 
l'imagette consid6r6e, Une valeur inferieure & ce seuil 
5 donne une portion insignif iante, telle que J3 . 
Cependant, les portions insignif iantes des image ttes ne 
sont pas reconstruites, c'est-S-dire que les calculs de 
rStroprojection ne sont pas effectues pour elles. 

On procSde concretement en construisant un 

10 masque numerique avant de retroprojeter l'imagette. Sur 
les details horizontaux Dh, le masque est construit 
pour la premiere fois pour l'imagette Dh3 . II prend une 
valeur egale k 0 pour les portions insignif iantes 
telles que J3 et Sgale a 1 ailleurs. Les coefficients 

15 du masque se suivent dans un ordre determine, par 
exemple ligne aprds ligne. Les calculs de 
r^troprojection sont appliques sur l'imagette 
consid6r6e dans 1' ordre des coefficients du masque. 
Quand un coefficient est egal & 0, aucun calcul n'est 

20 effectu6 pour le point correspondant de l'imagette, 
auquel on attribue une valeur nulle dans l'imagette 
retropro j ectee . 

En abordant 1' image suivante des details 
horizontaux (Dh2) , le masque numSrique permet de ne pas 

25 considerer la portion J2, dont les points ne sont pas 
trait^s par 1' unite de calcul qui effectue la 
r^troprojection. Si par exemple le masque de l'imagette 
Dh3 possede un coefficient nul a la ligne i et la 
colonne j, on pr^voit que les quatre points des lignes 

30 2i et 2i+l, et des colonnes 2j et 2j + l de l'imagette 
Dh2 auront aussi des valeurs insignif iantes . Les 
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calculs de retroprojection ne seront pas effectues pour 
ces points. 

Un masque numerique est ainsi construit 
pour chaque niveau de decomposition. Le masque 

5 decrivant l'imagette Dh2 comprendra par convention des 
coefficients egaux a 0 pour toute portion homologue k 
une portion aux coefficients nuls (comme J3) du masque 
de l'imagette correspondante a plus petite echelle ; 
pour determiner les valeurs 1 ou 0 des autres points du 

10 masque de Dh2 , on utilisera encore les comparaisons 
avec le seuil convent i onne 1 . D' autres portions 
insignif iantes, ayant des valeurs nulles de masque, 
pourront apparaitre. C'est ce qu'on a illustre sur les 
imagettes Dvl, Dv2 et Dv3, des decompositions 

15 successives des details verticaux. L'imagette Dv3 
comprend une portion insignif iante K3 dont on retrouve 
les homologues K2 et Kl sur les decompositions & plus 
grande echelle. L'imagette Dv3 ne comprend, dans cet 
exemple, aucune autre portion insignif iante, mais il a 

20 ete possible d'en trouver trois autres, notees L2 , M2 
et N2 , sur 1 ' imagette suivante Dv2 . Au niveau de 
decomposition a plus grande echelle, celui de 
l'imagette Dvl, on ne s'occupera pas de retroprojecter 
les points situes aux portions LI, Ml et Nl homologues 

25 de L2, M2 et N2 . 

On va maintenant examiner d' autres 
particularites du procede de 1 1 invention favorables, a 
une acceleration des calculs. 

La premiere repose sur l'egalite entre la 

30 transformee de Fourier d'une projection de 1 1 image Sl un 
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angle fix£ et la transf ormee de Fourier de 1' image sur 
une droite de meme angle passant par l'origine. 

En se reportant a la figure 3, ou un jeu de 
projections & ete convert i dans le domaine de Fourier 
5 pour dormer des projections de nature f rSquentielle 
dans le systdme d'axes notes £l et £2, le support du 
jeu de projection dans le domaine de Fourier comprend 
les valeurs comprises entre -v 0 et v 0 pour £ x comme pour 
£ 2 . Une decomposition en ondelette, comme celle de la 

10 figure 1, fait apparaitre 1 ' approximation encore not^e 
AA dans les frequences les plus basses, autour de 
l'origine £l = £2 = 0, alors que les details se 
trouvent de part et d' autre de cette approximation, £ 
des frequences de plus en plus grandes pour les details 

15 a grande £chelle. 

On a represents une droite passant par 
l'origine, en faisant un angle 9 avec l'axe horizontal 
£l . Cette droite traverse 1 ' approximation AA, ainsi que 
les imagettes des details verticaux Dvl et Dv2 . Mais on 

2 0 est assure que les projections faisant cet angle 0 

seront tout £ fait inutiles pour les calculs de 
rStroproj ection des details diagonaux Dd et horizontaux 
Dh puisque la droite d' angle 6 passe k 1 1 ecart de leurs 
imagettes. L' application de 1' invention comprend alors 
25 une selection, pour chacune des categories d' imagettes, 
des angles de projections qui seront utilises dans les 
calculs. Le calcul sur le support des projections est 
elementaire . 

La figure 4 reprend la division d'une image 

3 0 decompos6e en ondelette et transposSe dans le domaine 
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de Fourier. Une reconstruction parfaite peut Stre 
obtenue avec l'emploi d'un nombre determine de 
projections, dependant de la frequence (discretisation) 
de reconstruction de 1 1 image de l'objet. Ainsi, pour 
5 reconstruire 1 ' approximation AA, on peut demontrer 
qu'il suffira de choisir un nombre de projections, 
parmi celles qui ont ete faites, correspondant a une 
frequence maximale VI servant de rayon a un cercle 
circonscrivant la representation f r^quentielle de 

10 1' approximation AA dans la decomposition. De mime, les 
details a petite echelle seront compldtement rendus en 
utilisant un nombre de projections correspondant & la 
frequence v2 dans le cercle ayant ce rayon, dans le 
systeme d' axes £l, £2, et circonscrivant les 

15 representations f requentielles du groupe de details 
Dv2, Dd2 et Dh2 a meme echelle. Le support des 
projections dans le domaine de Fourier permet done de 
connaitre aisement les frequences maximales necessaires 
pour la retroprojection parfaite des imagettes 

20 respectives. 

La mise en oeuvre de 1' invention pour 
exploiter ces deux particularites comprendra done, a la 
retroprojection de chacune des imagettes, une selection 
des projections utiles pour cette retroprojection, en 
25 ecartant les autres ; et eventuellement une restriction 
du nombre de projections utiles ef f ectivement 
exploitees par le calcul pour n'en garder que le nombre 
utile . 

La figure 5 est un organigramme de 
30 1' ensemble du procede decrit ici. 
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RE VEND I CATIONS 

1. Proc6de de reconstruction d*une image & partir de 
jeux de projections de cette image, comprenant 
5 successivement : 

- une s6rie de decompositions successives 
des jeux de projections en ondelette dormant des 
imagettes des jeux des projections, comprenant des 
images d f une approximation (AA) et des series 

10 successives (Dd, Dh, Dv) de details homologues de 
chaque jeu, 

- dans chacune des series et successivement 
pour des imagettes de details ayant un nombre croissant 
de points, une recherche de portions insignif iantes 

15 jugees pour depourvues de contenu, et une recherche de 
portions homothetiques des portions insignif iantes dans 
chacune des imagettes de details qui suivent dans la 
serie ; 

des r^troprojections des imagettes des 
2 0 jeux de protection en imagettes de 1 1 image a 
reconstruire , avec une omission des r§troprojections 
pour toutes les portions insignif iantes et toutes les 
portions homo t h§ ti que s, 

et une combinaison des imagettes k 
25 reconstruire par une inversion de decomposition dormant 
ladite image . 

2. Procede de reconstruction d'une image suivant la 
revendication 1, caracterise en ce qu'il comprend une 
selection de regions d» angles (0) des jeux de 
30 projections qui sont utilises dans les 

rStroprojections . 
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3. ProcSde de reconstruction d'une image suivant la 
revendication 2, caracterisg en ce qu'il comprend une 
selection d'un nombre de projections qui sont utilis^es 
dans la ret roproject ions . 
5 4. Proc6d6 de reconstruction d f une image suivant la 
revendication 2, caracterisS en ce que la selection des 
regions d' angle est r^alisee en fonction d'un support 
d f une transform^ de Fourier des jeux de projections 
decomposes en ondelette. 
10 5* Proc£d6 de reconstruction d'une image suivant la 
revendication 3, caracterise en ce que la selection du 
nombre de projections est realis£e en fonction de 
frequences maximales d'un support d'une transform^ de 
Fourier des jeux de projection decomposes en ondelette. 

15 
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FIG. 2 
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MESURES (JEUX BIBIMENSIONNELS 
DE PROJECTIONS, OU PREMIERE 
IMAGE DE L'OBJET) 



REAMENAGEMENT DES 
MESURES CONIQUES EN 
MESURES EN EVENTAILS- 
PARALLELES OU DES 
MESURES EN EVENTAIL 
EN MESURES EN PARALLELES 



DECOMPOSITION EN ONDELETTE 
DE L' IMAGE EN IMAGETTES 



SELECTION DES 
PROJECTIONS UTILES 
POUR CHAQUE IMAGETTE 



FIG. 5 



IMAGETTES DES DETAILS A 
LA PLUS PETITE ECHELLE 



DETERMINATION D'UN 
MASQUE DES IMAGETTES 



RETROPROJECTION 
DES IMAGETTES DE 
MEME ECHELLE 



PASSAGE A 
L'ECHELLE 
SUPERIEURE 




NON 



SYNTHESE DE L'IMAGE 
DE L'OBJET PAR COMBINAISON 
DES IMAGES RETROPROJETEES 
EN ONDELETTE INVERSE 
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FIG. 7 




